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I. EQUATIONS OF SURFACE SHIP MOTION 


A moving ship is a body with vix degrees of freedom. 
These degrees of freedom are generally chosen as follows: 
a. Linear displacements along the three axes 
through the center of gravity. 
a.l Surge -- along X axes 
a.2 Sway -- along Y axes 
a.3 Heave ~- along Z axes 
b. Rotations around the three axes through the 
center of gravity. 
b.1 Roll -- around X axes 
b.2 Pitch -- around Y axes 
b.3 Yaw “~~ around Z axes 
Further reduction in the complex nature of the equa- 
tions can be brought about by choosing an orthogonal 
axis system parallel to the principal axes of inertia so 
as to eliminate products of inertia in the motion equa- 
tions. For practically all ocean vehicles, with extreme- 
ly few exceptions, a longitudinal axis (X axis) in the 
centerline plane, a downward (toward keel) axis (2 axis) 


perpendicular to the X axis in the centerline plane, 
and a 


as ‘ 


centerline plane satisfies this requirement. 


ee re EN Re SE ES ; a es PE ya CLC RS “S35, =F: wee a 
FS RES Re TI, ETE RAE ear Ay IER LEY PERE TEA er PIR RET RB, SEO REST RE RE CT FEE TER RE ART POP BT NEE PET NIE BT SBE RE a 
: mF i 


FERS meron <n ki eis’ 2x “ed 


~HEAVE 


=- SWAY 


SURGE + 


+ PITCH 


+yYAw ’ 
+ SWAY 


+ HEAVE 
Figure 1. Surface ship in six degrees of freedom. 


For the exceptional vehicle which has a very peculiar 

and significantly large asymmetrical mass distribution, 
; it is necessary to include the products of inertia. 

The X, Y%, 2 axes form an orthogonal right hand system 
of axes fixed in the vehicle. The axes and the associated 
components of the pertinent physical quantities are 
defined below: 

The longitudinal X axis (in the plane of symmetry) is 
positive in the forware direction, usually parallel to 
the keel or calm water line. I£ the upper and lower halves 
of the body are symmetrical, then the axis is the inter- 
section of the two planes of symmetry. 

The Y¥ axis is the transverse axis perpendicular to the 
piane of symmetry and positive to the starboard. 

The Z axis or downward axis in the plane of symmetry 


(X,2) is perpendicular to the X axis and positive downward 


towards the keel, 


64 Ry OTR SA A 


SERIE RTS ae wea re eign yn et pee rer PN ERNE EN TIO eH POET ARMREST A Tem HY MS SHS ABT 


ce 


AMIE GIG ARS By A gt oie 


BT FM Ta, 


FE 
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EAT ERT STE ALY GU TR TTT TA ETT OTT MIE AE sae 
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email xine 


o 


dha 
fe 


nan A A 


i,j,k unit vectors along the X,Y,2 axis respectively. 
R xX,y,% vector distance of a point from the origin O, 
and the corresponding components along the X,Y and Z axes. 
R = ix + Syg + keg 
u u,v,w velocity or the origin O (on the body) and 
the corresponding components along the X,Y and Z axis. 
t= iu+ jv + kw 
$ p,q.xr angular velocity of the body about the origin 
and the corresponding components about the axes. 
$= ip + jn + kr 
The moments of inertia of the body about the X,Y and 


Z axes respectively I, Ly I: 


y 

F X,¥,Z, force acting on the body and the correspond-~ 
ing components along the axes. 

f= ix + jy + kz 

m K,M,N Moments acting about the axes. 
hx GK + 5M + <N 

Newton's law of motion for a rigid body can be written 
as two equations, one a force equation and the second a 
moment equation provided an origin is taken at the center 


of gravity and the axis system is fixed in space. The 


equations are 
> 


= _d J a > 

PS ae (momentum) ar (m Ug) 

mad 1 aaj S268 
m= a5 (angular momentum) = ar (2 ) 


eS eae ee eee EES WEI eee ta rare oe ee oe nA ry, 3 a 
ipa ch aah aS : FE SEE EF NEES HURTS SN EE ELE FOTOS OR WELL TE OE y OR EE EIT PES EPL BEE NS OE Ie EOL LIEE LIE TIS MRS TATE TTS OLE LTO TG 
SE: s i SE PLE PS TE ETT OTE AT OL TT ETT 


— wm 


where the subscript G refers to an origin at the center 


of gravity and m is the mass of the body. For a mass 


MPT 


essentially constant in time 


For an origin not at the center of gravity of the 


UE AE AS ORE RTO TE PTT Nh 1 


body and in a system of axes fixed in and moving with 
_the vehicle. 

> > 

q = 

Ug =u, +h xR, 
where UL is the velocity of the origin in space. However, 
since the origin is on the surface of the earth and the 


earth rotates, then 


> > > 
BU, =U+& x R, 


where f is the geographical velocity of the body, om is 
the angular velocity of the earth, and R, is the radius 
vector from earth's center to the vehicle. The force 


equation becomes: 


> d , > > 
F=m ae (U + ae xR +2 x Ra) 


a +@x R,) +mid_x Rt8 x R] 


b>) a 
ui 
3 


=m (8 +8 x R,) +m (Sx rs +0 ox x R,) 


; + > 
since % = 0 and Ry, a eR 


| 
C+ 
+ 
Ot 


> 
the term mi x is the coriolis force and mi x Cm x Ry 


is the centripetal acceleration due to rotation of the 


rag 
e PAM ae mene 
a 


earth, These two terms are negligibly small when com- 


pared with the other forces, then 


F =m oti +x Rg) 


Finding the derivatives of unit vectors (change in 


direction) 
where di = -kd@ di = ja¥ di = 0 
dj = 0 aj = -ia¥ aj = kab 
dk = ia 8 ak = 0 dk = -jab 
Adding the contributions 
di fy 4 34H. fae 
ae" 2°9 + Jae ~ Xa 
d od . “d 
at = ~ige + 3:0 + KEE 
ak - dé <a - 
ae = i ge gg + k+0 
Q = ip + iq +k 
= ot = 38 a Wy 
Pr ae ' 9" ae * * ae 
oe i jk 
= y i= > gq r 
dt 1 oO O 
t= iu + jv +k 
> : 2 - 
Ro = iX,+ jY¥Q+ K2a 
* A . " 2 a 4 
a = id + uj Se + jv + v2 + kw + woe 
= ie $ (add = 7.08, its 284 et soy as 100. 
ae Wae. mate. 7 1 *dt dt 
es *d8 = 4d 
+ kw + wise Jae) 


= iut u(jr- kq)+ jv+ v(-irt+ kp)+kw + wliq - jp) 
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2.2 Kk 
a x Ro = ep q-° = i(qZ,- r¥,)- j (p2,- rX,) + k (pY,~qX,,) 
*e 4g. 46 
Se (ixB) i(Z.q-¥or)+ j (Zgar-¥or?)+ k (Ygqr-2,97} 
e ° ° 2 2 - 2 
j(X, r- ZqP)t+ k (X gtP ~ Zap )+ i(Z,rp ~ Xor ) 
2 fs 2, 
k (gb XQa)+ i (Y¥gps qq )+ 3 (Xepd ~ Yop 
-  _du.a + 
from F = ix + sy + kz and F = mae + ae x Ro) 
Rewriting and grouping all i terms equal to X, all j 
terms equal to Y and all k terms equal to Z yield 
X = m (u-vrtw +Z.q + rpZ -r7x + pqY, - 2y ) 
9t4,9 P4, Gg Pata ~ 246 


m(itwg-vr-X, (r7+q°) + ¥, (pq-r) + 2, (rp+d) 


(1) 


Y= m(vtur-wp+X, (pqtr) - Yq (pete?) + 2, (pq-x) 
Z = m(wivp-uqiX,(rp-g) * ¥,(qxtp) - Z,( 24 q”) 
& mMi(wrvp Girag iEP Gg GstenP. 4g‘'P 7 
F m, = S-(angul tum) , 
rom M, = g¢langular momentum 
d 
= Sci, pt sy q+ kt, r) 
Xg 26 


G indicates an origin at the center of gravity 


a 7 2 2 

ae I, ™ (¥. + Ze ) 

x _ 2 2 

we Iy m™m (2, + Xo ) 

rps 2 ee ee a 

Zo z G G 
> > > > +> > > > 
M M, + Ro x F or M, = M - Ra x F 
11 
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After manipulating and using the results for the 
: derivatives of unit vectors (same as ahove) expressions 
for K, H and N are obtained. 


K 


Typ +(I,- Ty)artm[¥,, (w+ PV - qu)- Z6(v+ru-pw) ] 


er ree Reet 
St ERE Leeann an ev oe ay SRE AME AeA ae 


= Iyq +(I,- I) rptm[Z,, (ut qw ~ rv)- X, (wtpv-qu) J (2) 


= 
I 


TMC CL ae Te HUTS 


x +(T,- Ty) pqtm[X, (vt ru - pw)- YX, (utqw-rv) ] 


z 
i] 


Ty 


The terms (qw-rv), (ru-pw) and (pv-~qu) are gyroscopic 
effects. ; 


The relationship for forces and moments can be ex- 


pressed 
Force) 
) = £(properties of body, properties of 
Moment) motion, properties of fluid) 


SOT LT EI ETN EET LT EE NTT eH TT a 


2 ae F ree er 
SORT TE ELAR SRC eH 
. Ss bed ps OAS whe ak em S.A RNY PURSE Sree ie eM SEE coeURCon Co eee aye we Rewer A me 


For a particular ship, in a given fluid with no 


excitation force - so 


f£(properties of motion) 


Ry My 
ee ee Se 
tt 


E(X 1X 5050018 We UrViWePeeL 78, 
Prdet 5,5", 8") 
The Taylor series which has the following form may 
now be applied to linearize the equations about an 
operating voint Xo 


C(x) = 2K) + KS £00) 4 (x-R)? a 2 4 
& Fe Be Nee Se * “oO ax Ay oO sere S Ae ihe ee 


Appi, this to EUX V2) we siwees 
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For the case of £ (properties of motion) let (X-%) 


= AXSs (¥-¥,) = A¥, and (2-2,) = A, 


terms second order and higher are neglected for small 
perturbations. 


From X «qguation the linear terms are obtained: 
= ; bf : 
Caos t ( = dat AX (az ee AY, CF )+ AZ (Fz, on oe eAV v (pe) + vee 


Tr. defining relations are: 


(é f) = (3 =) =X 

suo Sul uau, u 
eH WweWo =X 
Aw = (ww) =W, bed 0 


Au = {u-u,) 


The force equations then become 


X= X, + Xy %, + X Yi eee oX,0 + coe KO, 


; Xo re) 6 

5 a oie enc Yy Xo + YY +244 .¥Q8 + reer YA (3) 
oO °o 

Zz = Zo + ey Xo + ae Yo $444 6Z90 + vee eZ dy 


° 
A similar derivation can be done for the K,M and N equa- 
tions. The preceding X,Y and 2 equations may now be 
equated to the linearized X,¥ and Z (equations (1)), e.g., 
fox the Y eguation without roli, pitch, and the center 

of gravity at Xo= 0, Yor 0 and 2g 0 gives the linearized 
equation. 


Yr m(v + ur) 
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+ “1 ey oye r 
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Xo 
Expressions for X,2 and K,M and N cas be determined in 


a similar procedure. 


In order to obtain the equations in a non-dimensional 


form some definitions will be given, and applied to the Y 


} 
& 
é 
7 
t 
F 
: 
E 
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E 
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i 
h 
F 
f 
a 
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force equatious as an example of the process. 


Froude number = we 
vg 
u'.v',w' = pa 
¥g/% 
ts t(v¥g/2) 
X'oy',Z' = X,Y ,2% 
, ¢ 3 
pgk 
Cet er X,Y,2% 
z ou g 


After replacing and adding the eiiect of waves, the 


Y equation becomes: 


vit rlal « 5 ne ay: + Cy P! + Cy a Cc ér) 
v p x 


+ (YoeP! te Yoer' + Xi) 
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€ II. DIGITAL COMPUTER SIMULATION 


The six equations of motion after rearranging by plac- 


va BD 


nah 


ing the second order terms to the left and the rest on the 


right side becone: 


awit arin, 


aaAt+baBt+caCtdaD+. aE+faF = ~ (a a, Ata a,Atb, a, Btb,a,B 


+ C,a,Cte,a,C+d, a, Dtd,aoD 


+ €,8, Ete ,a,Etf a F+f,a5F) 


+ IF1 
abA+bbB+chC+dbD+ebE+£bF = ~ (a,b, Ata ,b Atb,b, B+boboB 


+ cb, Cte gPoC+d, by Dtd bod 


FFE gi TOTP TL TEE ONE TENT 


* 


bl oe Si 
a3 TS OP PRES AI EO GRRE ERIN CTI 27 


+ eb, Ete bo Etf b  F+f boF) 


+ IF2 


acA+bcB+ccC+dcDteck+feF = ~ (a,c, Ata,c,Atb,c) Btby cB 
+ €,)¢,Cte,coCtd,c,D+docD 


D+F 
2egEt£ co P+f clF) 


+ IF3 


adA+bdB+cdC+ddD+edE+faF = ~(a,d,Atad Ath, d, B+b,dB 


+ 67d, C+e,d,C+d, dy 


+ e,d, Ete,d Erfjd)F+f4d, 


+ IF4 


D+é,4,D 


F) 


aeAt+tbeB+ceC+deDt+eeE+feF = ~(aje,Ata,e,Atb) 6, Btb,e 5B 
+ ce, Cte ,e,Ctd,e,Dtd,e5D | 
+ e)e) Ete,e,Ftf eu Ftf,e5F) 


+ IF5 


ee anne te 
Mat, Bane Soy sn sae oy 2 . 
WUE te tet Fein ee an a ee 


a Sa ai seo Hh AE See aS Fe ER ae a Le ene ee en eer ee SEA CREST TERRIAR CREME PUTA RUE OT oe ONO LEB Sy epee FT ES se Et a OR NE | Guear > 


eae ze 


afAtb£BtcfC+d£D+efE+fEF = ~(a,f£,Ata,f,At+b, £)Btb, £28 


+ ce, f,Ctc,f,Ctd,f,D+d,f,D 


+ e,f Ete f Ett f Ft+f,f5F) 


+ IF6 


where A=u, A=u, B=¥, B=v, C=w, C=w, D=p, D=p 
E=q, E=q, F=r, F=r, terms IF include all non-linear 
terms such as wave force, wind, rudder deflection---etc. 


In the six equations, the right can be set equal to 


I: I, atta I respectively, thus: 


an 
H 


= ~(aja Ata F) + IF1 


2 
2 = ~fayby Ata boAtby by Bt ---<-----~ fob. 


a,Atb,a, Bt Stet teteate tes fia 


+ 
! 


FE) + IF2 


It 


~ (ajc) Ata,c a+b, c,Bt ner ment aa f 


~ (ad, Ata ,d,A+b,d, B+ Ao a a te em f,4,F) + IF4 


F) + IF3 


aad 
H 


5 = ~ (aje,A+a,eoAtb) e, Bt SSSSssea= “fi,e,F) + IFS 


I, = ~(a,f Ata f.,At+b f.F) + IFG 


7 ay ahd) Nes aaa aa f 
the equations thcn have the torm that follows, 
aaA + baB + caC + daD + eak + faF = I 
abA + bbB + chC + dbD + ebE + fbF = I 
acA + beB + ccc + deD + ecE + foF = I 
adA + bdB + cdC + ddd + edE + fdF = I 


aeA + beB cec + deb + eeE + fer = T 


+ 


afA + b£EB + c£C + G€D + efE + FEF = 1 


3 ROE REE FE ; 
ae % F Qe See TEE ay 


Sahm en fee Sen ee ce eee oe 


bi STE ERE TEL SRF Aa” flea 
. ee fiat 


expressing in matrix form: 


aa ba ca da ea fa A I, 

ab bb cb db eb fb B I, 

ac be cc dc ec fc Cc I, 

ad bd cd dd ed fa D| “| x, 

ae be ce de ee fe E I. 

af bf cf ct ef Ef F I¢ 
Apply Cramer's rule to solve for A, Bo-wnn- eee F in 


terms of Ty cael Te 


Ty ba ca da ea fa 
I, bb cb db eb £b 
I, be cc dc ec fc 
I, bd cd ad ed fa 
I, be ce de ee fe 
I, bf cf df ef ££ 


eet OES tt SNR, 


aa ha ca da ea fa 


A= 

ab bb cb db eb fb 

ac be cc dc ec fc 

ad bd cd dd ed fd 

| se be ce de ee fe 

| at bf cf daft ef ff 
define the denominetor determinant : 4 and for the 
A 


oo, £ am 
GUne wu 


nominator let cofactor of ty 


A cof. ab, -~ere---- and coefactor of I 


cofactor of I, 
é cof. af equations becoines: 


(cof.aaI,+cof.abl,+cof.acl,t+cof.adI,+cof.ael,+cuf.afl,) 


A = SS ee 


_ a <a 


4 
4 w 
he 


2S Ee ag” Meester Tne ng mat RCE 


RIMES Bung Ho ARR a SETS RP 


RE Sf AD MI NTR 2 OPE PLO ENT IM! EFT MD TOI Es PLATT 4 2, 
- 


POE LSE TRE PETE FIL YT TORY 
sf 


a = a = aoe es 
ES TE OE CREME SS CEL ee eT TY Ea, ae OS EPR SRILA AY aa Te ee ae 2 
PEER : aS as Ne eS ee re ae 


Ir. the same way solve for B, C, ---------~ 


oe 


{cof bal, cof boi, cof bel, cof bal, cof bets cof bet? 
A 


(cof cal, cof ebi, cot col, cof cdl, cof cel, cof cfI,) 


c= — 


(cof dal, cof dbI, cof del. cof dat, cof dei, cof dfI,) 


B= 

ae (cof eal, cof ebl. cof ecl, cof edl, cof eel, cof efI,) 
A 

oe (cof fal, cof fbi, cof fel, cof Fal, cof fel, cof ffI,) 
A 

Then the value of A, A, B, B --------~ F, F by integration 

such that 


A block diagram to compute ail of the variables in the 
set of equations is presented in Fiq. 2. 

In the computer program that is used for simulation 
a.\l six equations for six degrees of freedom are used, but 
are interested in less than six degrees of freedom. The 
same program can be used by setting the coupling terms of 
non-used equations equal to zero and one in terms of princi- 
pe agonal, 2.4. only three degrees of freedom are used 
i 


n this study, surge, sway and yaw, then all coupling terms 
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Fig. 2 
TRANSFER FUNCTION BLOCK DIAGRAM SIX DEGREES OF FREEDON(EQn.A) 
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TRANGEER FUNCTICO.; BLOCK DIAGRAM SIX DEGREES OF FREEDOM(Eqn.B) 
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Fig. 2 
TRANSFER FUNCTION BLOCK DIAGRAM SIX DEGREES OF FREEDON(Eqn.C) 


TRANSFER FUNCTION BLOCK DIAGRAM SIX DEGREE OF FREEDOM(T@n.D) 
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TRANSFER FUNCTION BLOCK DIAGRAN SIX DEGREES OF FPREEDCw(Eqn.£) 
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Se Bn hd pest a alee 
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BFE RE ES” ARTE RES PT aS DT EOP EN AD UIE PND HE. e. sas F we 
; PEPE S BO EN A ILI ERE GALT TET OT ERE EEE TE OT ROTA SOE TLE FLEES ST TT FTE Beam: 
AE cee es START FTE OYA SL BAS ES Ee CET 
3 iy 


TAB, 8 
re te 


we yd eh wD TI 
eg) 
arth Neal hoe 


eed 


oa 


are set equal to zero and unused terms on the principal 


Sea ore eee 


; diagonal equal to one, for example: 

Ee 

aa ba 0 0 0 fa A T 

i ; ab bb 0 0 O _ fb : i, 

; 0 0 1. <G: @ 0 : | 13 

Fe 0 0 0 1 0 0 D ly 

0 0 0 0 i: 0 E Te 
at bf 0 0 OO f€£ | :| Le 


and the left side of the unused equations are set equal 


to zero. 


With this program non-linear terms can be added such 


as, rudder deflection of waves and wind, etc., which will 


be done by adding all of these whose sum is IFy +9. 


LE Bea ett Ste een otarree eper ance enartyon 


1F) = KA,x Dr+KA.x Ds+KA,x Db+NA 


where Dr, Ds and Db are rudder deflection, canard de- 
flection.......... ete. NA is the sum of all non~iinear 
terms that effect the surge equation (X equation). 

The program that will be used for solving these equa- 
tions is the "Continuous Systems Modeling Program" (CSMP) 
{Ref. 3] in which all constants are declared in the first 
section and then set the value of matrix for aa, ab, ac 
seeseee and so on (in program AAA is used for aa, AAB for 


“ 


mn the int 1 eaction values of 


ab eoeeeean ALF Lox Iij. In set 
the COFACTORS aa, ba....are determined. All of the CO- 
FACTORS and the subprogram VALUE is used to compute. This 


subprogram finds the value of the determinant of the 
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Fe ee ne ee OER ET EERE TN ae tee gee Maa Pe ay S| 
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ee . 


SOE TS PA EY POTEET 


matrix. For all of the COFACTOR terms the element is set 


equal to one and the rest of the elements in that row and 
column are set equal to zero. For example, to find the 
value of COF.aa the following array is obtained: 

1 0 oO oO O 9 

0 bb ch db eb fb 


COF aa = be cc ae ec fe 
bad cd dd ed fd 


be ce de ee fe 


> i > ee = 


bf st df ef ff 


LL AE RR EN RONEN LT LO Te ON RINE NE AA LORE HORT AY 


(In the computer program BAP. is used for a,ay, GAA for 
AoAoeee)- After the value of A and all cofactors are 


determined, the dynamic section is used to determine BAA, 

RAB ......GAh, GAR ._...., (if those terms contain variables). 
In the dynamic section all variables that are functions 

of time are determined. The defining relations of the 

variables are also included in the dynamics section, i.e. 

UDOT = ADOT (U=A), U=ADOT.......etc. XH,YH,ZH are deter- 


mined and are the vector terms, X,Y,Z whose origin is ‘fixed 


on the earth (relative to the earth). 


Pe. mr et, ee ae me er TE me a a a a te 
CREPE EITN IY IESE NRE TEE RN NOC EMER RARN ENCE” NNT RRR HMM Ie Em eemnawertNs Weeu SteNnie ERENT = areca e 
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III. STUDY OF SHIP "D" PERFORMANCE 


In this section the computer will be used to solve 


the equation of motion describing ship "D",. The hydro- 
dynamic coefficients and constants that were obtained 
from NSRDC (NAVAL SHIP RESEARCH AND DEVELOPMENT CENTER) 
{Ref. 6) fur this study will concern only three degrees 
of freedom such that SURGE, SWAY, and YAW. 

Equations of motion of the ship (nondimensional) 


Axial Force 


m(u-vrtwg) = 5 £(Xo09° + yor + XptP) 
+ 5 (x30 + X vr + Xygvd) 
+ Brix") 
+ FeO (gp gy St HX 5 55 558 *Kgn gp OP) 
+ xr (au? + Aureus + Au”) 
Lateral Force 
m(vtur-wp) = 5 L(Y sp + Yor + Yoghd) 
+ 5 (YpWPt¥y | p]v ie] +ay r+ jv+uy p) 
+ 5 (u¥v + Yow + Yiy}ylvly) 
+ ort u*(¥, 6x) 
Yawiny Moment 
Ir + (Iy-Ty)pq = 5 (Nex + Nesp + NygPd) 
+ 5p (N,VEUN, pruN THN, WPAN) 1 |v] x) 
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Fe TEMG are wees oe fee me eee 


wr 


9 
+ oe aa ls |v}v) 


Iviv 


fe) 2 
+ —— u” (N, 6r) 
; 227 6x 


p, the density of fluid is taken as 2 and the terms 
including w, p. q (heave, roll and pitch) are set equal 
to zero. 

The nonlinear terms such as the squared terms and 
product terms of v and r are omitted initially. This is 
in agreement with the small perturbation theory. 


After rearranging, the equations become 


2 
OF 2u 
° . a 0" 
(Ys - mv + Yur aaa: YVv - uy or = ae Y,,6r 
* 2 
Nyy + (Nye-Ip)x = = NV- SNC - oa Ng Sr 
2 2 R 


Set the left side of the equations equal to Tjr---I¢, 


then the matrix equation becomes 


(Xm) 0 0 oO 0 0 i) 
0 (Y.-m) 0 0 oY, B I 
0 0 2 0 0 0 C Is 
0 . 0 0 1 0 0 D Iy 
0 0 0 0 1 0 E I. 
0 Ns/l 0 0 0 (Ns-I,) F I, 


where 0 =A,v=B, r=F. 
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Also set the right side of the equations equal to I, 
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All Zero IF4 


where u =A, Vv = B, x = F, and IFI1 = 


nN 


= = eh 
IF3 = 0 
IF4 = 0 


IF5 = 0 


Z 
u 
rot ey 


IF6 = -N 
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A. NONLINEAR TERMS NOT INCLUDED 


TABLE I 


Hydrodynamic Coefficients and Constants ot Ship "D" 
for Linear Terms ‘Non-dimensional) 


m = 09,0045 


u 


Iz 0.0003 


Nr Q.0012 


i] 


Nr = -0.0002 
Nv = -0.0012 
Nv = -€.0601 
Xu = -0.00036 
Yv = -0.0025 
Yvy = -0.,0063 
Yr = 0.004 

Xe, = 0.0011 

XY = 0.0019 

Ne = -0.00084 

és anc 6b equal zero 
In the computer prosgram, set all coefficients in 
secticn 1 and set AAA = Xa-m, AAB = 0 ---- AFF = Nr-Iz. 


After this, use subprogram value to find the determinant 


and coefficient of AA, BB --~- and then set BBB = 7 Yv , 
BFB = uYr ------~ BFF = z NE : 
2 
IFl = KAlér where KAl = -X,_,. 6x aaa 
IF2 = KBlér where KB2 = oY u7/L 
2,,2 


IF3 = KFlér where KF1 = -N . 
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Following is the program that used CSMP to determine 


the turning circles for rudder angles of 15° (~0.2619 rad.), 


Ce tein eo a ar re in ee 


25°(-0.4365 rad.) and 35°(-0.6111 rad.). Result of this 


PA TOS SE I ARR a ay 7! 


cos 


study is presented in Fig. 3. The turning rate as a func- 


tion of time is interested and the resulis of this analysis 


b 
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are presented in Fig. 4. The ships turning performance 
expressed in transfer ship lengths as time is shown in 
Fig. 5 and the heading angle as a function of time is 
given in Fig. 6. Fig. 7 shows the results of the zig-zag 
maneuver, curve shown the yaw angle and rudder angl2 in 
degree as a function of time, for this study the same pro- 


gram was used, but set DR in dynamic section: 


, 7 
EUTELSEA 2 2S givers > ee peter tas wy a Baer ste YET aN One rah 


~RAMP (65.0) ~RAMP (145.0) +RAMP (170.0) +RAMP (250.0) 260. 
~RAMP (265.0) -RAMP (345.0) +RAMP (425.0)+RAMP (440.0) ) 


and use prepare statement prepare X, YAWD and prepare X, 


DOO (YAWD = YAW* 57,273, DOO = DR* 57,273). 
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The computer program 2 is the same as the computer 
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program 1, but nonlinear terms are added by setting: 


ena 


NA = NAl + NA2 + NA3 + NA4 
és _ 2 2 
where NAl = 1(4,49 + Xe + Xp pXP) 
f NA2 = ~(mvx + X vr + Xt + mwq) 
oe 2 
NA3 = (Xy¥ )/i 
NA4 = -(Alu? 


2 
+ A2uu, + A3u, fl 
NB = NB1 + NB2 + NB3 


where NBl = -1Y 
pg?4 


; NB2 = -(¥,, 5¥P + Yyjrpviel+ mur + mwp) 
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NF = NF1 + NF2 + NF3 


where NF1 = “Ni gPa + (Iy-Ix) pq 
NF2 = - (N,,54P + Niyjrlviz)/2 
an 2 
NF3 = -(N\ wv + Niyjylvlv)/2 


Again setting terms that include w, p, and q (heave, roll 


and pitch) equal to zero, and set number of known terms 
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TABLE II 


Hydrodynamic Coefficients of Ship "D" for Nonlinear 
Terms (non-dimensional} 


Kor 0.00005 


i} 


0.00241 
~0.00341 


Ld 
1 


VV 


Yeiv|> -0.0416 


Nyiv|= -&.01002 


PROPULSION RATIO 4 n 


n > 0.45 -1.0 < n < 0.45 n< 1.0 
Al ~0.00004 -0.00032 -0.00117 
A2 -0.00035 0.00070 -0.00100 
A3 0.00099 0 -0.00085 


Fig. 8 - Fig. 12 are the same results as Fig. 3 - Fig. 7. 
The results of computer program 2 are more accurate than 
the computer program 1, when compared with free running 
model test of NSRDC (Ref. 6]. Fig. 13 and Fig. 14 when 
studying the stability of the ship by applying a force 
moment to the ship in computer program, set 

NFl 0.0001* (STEP (10.0)-STEP(10.01)}) (can use NF1l be- 
1 ause NFl in this program equal zero) 

Study direction stability of the ship by plotting 

direction of the ship (used advance VS. transfer ship 
length)and check heading angle of the ship by plotting 


YAW VS. TIME. 
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Iv. CONCLUSIONS 


The equation of motion of surface ship and computer 
program developed here including all of six degrees of 
freedom, but the study in III concerns only three degrees 
of freedom (surge, sway, yaw) because hydrodynamic coeffi- 
cients are not available; when the state of the art 
reaches Lhe stage in which hydrodynamics coefficients are 
available, this computer program can be used in all six 
degrees of freedom. 

Some results from III are not too perfect because the 


lack of some constants and coefficients such as the value 


of mass (m), initial velocity (ADOTO), command speed (UC), 
etc. But for study can adjust from curve for the model 
test [Ref. 6]. 

This computer program did not include some external 
effects such as effects of wave and wind, but these effects 
could be included in the program by adding terms to the IF 
equation. 

The following implementations are suggested for the 
future work. 

A. Study all six degrees of freedom. 

B. Study fox the effects of waves and wind. 

C. Study for control of the velocity and direction 


of the ship (by use of "MACROS', "PROCEDURE" or 


subprogram in CSMP). 
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